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Abstract. The structure of the F5 algorithm to compute Grobner bases makes it very efficient. 
However, it is not clear whether it terminates for all inputs, not even for "regular sequences". 

This paper has two major parts. In the first part, we describe in detail the difficulties related 
to a proof of termination. In the second part, we explore three variants that ensure termination. 
Two of these have appeared previously in dissertations, and ensure termination by checking for a 
Grobner basis using traditional criteria. The third variant, F5+, identifies a degree bound using 
a distinction between "necessary" and "redundant" critical pairs that follows from the analysis 
in the first part. Experimental evidence suggests this third approach is the most efficient of the 
three. 



1. Introduction 

The computation of a Grobner basis is a central step in the solution of many problems of 
computational algebra. First described in 1965 by Bruno Buchberger [7], researchers have proposed 
a number of important reformulations of his initial idea [5,6,8,9,15, 18 , 23J . Faugere's F5 Algorithm, 
published in 2002 [16], is in many cases the fastest, most efficient of these reformulations. Due to its 
powerful criteria, the algorithm computes very few zero-reductions, and if the input is a so-called 
"regular sequence", it never reduces a polynomial to zero (see Section [2] for basic definitions) . In 
general, reduction to zero is the primary bottleneck in the computation of a Grobner basis; moreover, 
many of the most interesting polynomial ideals are regular sequences. It is thus no surprise that 
F5 has succeeded at computing many Grobner bases that were previously intractable [I3|I16|. 

An open question surrounding the F5 algorithm regards termination. In a traditional algorithm 
to compute a Grobner basis, the proof of termination follows from the algorithm's ability to exploit 
the Noetherian property of polynomial rings: each polynomial added to the basis G expands the 
ideal generated by the leading monomials of G, and this can happen only a finite number of times. 
In F5, however, the same criteria that detect reduction to zero also lead the algorithm to add to G 
polynomials which do not expand the ideal of leading terms. We call these polynomials redundant. 
Thus, although the general belief is that F5 terminates at least for regular sequences, no proof 
of termination has yet appeared, not even if the inputs are a regular sequence (see Remark 1221) . 
On the other hand, at least one system of polynomials has been proposed as examples of non- 
termination (one in the source code accompanying [24 ]), but this system fails only on an incorrect 
implementation of F5. 

Is it possible to modify F5 so as to ensure termination? Since the problem of an infinite loop is 
due to the appearance of redundant polynomials, one might be tempted simply to discard them. 
Unfortunately, as we show in Section [31 this breaks the algorithm's correctness. Another approach 
is to supply, or compute, a degree bound, and to terminate once this degree is reached. Tight 
degree bounds are known for regular and "semi-regular" sequences [2j[20], but not in general, so for 
an arbitrary input it is more prudent to calculate a bound based on the data. To that end, 

l 



2 



CHRISTIAN EDER, JUSTIN GASH, AND JOHN PERRY 



• [T7] tests for zero-reductions of these redundant polynomials fSection l4.1lk whereas 

• pQ applies Buchberger's 1cm criterion (or "chain" criterion) on critical pairs (Section 14. 2[) . 

These approaches rely exclusively on traditional criteria that are extrinsic to the F5 algorithm, so 
they must interrupt the flow of the basic algorithm to perform a non-trivial computation, incurring 
an observable penalty to both time and memory. 

This paper shows that it is possible to guarantee termination by relying primarily on the criteria 
that are intrinsic to the F5 algorithm. After a review of the ideas and the terminology in Section [21 
we show precisely in Theorem [53] of Section [3] why one cannot merely discard the redundant poly- 
nomials in medio res: many of these redundant polynomials are "necessary" for the algorithm's 
correctness. Section 14.31 uses this analysis to describe a new approach that distinguishes between 
two types of critical pairs: those that generate polynomials necessary for the Grobner basis, and 
those that generate polynomials "only" needed for the correctness of F5. This distinction allows one 
to detect the point where all necessary data for the Grobner basis has been computed. We then 
show how to implement this approach in a manner that incurs virtually no penalty to performance 
f Section 14. 4p . Section 03] shows that this new variant, which we call F5+, 

• computes a reasonably accurate degree bound for a general input, 

• relies primarily (and, in most observed cases, only) on criteria intrinsic to F5, and 

• minimizes the penalty of computing a degree bound. 

Section [5] leaves the reader with a conjecture that, if true, could compute the degree bound even 
more precisely. 

We assume the reader to be familiar with [16] , as the modifications are described using the 
pseudo code and the notations stated there. 

2. Basics 

Sections I2.1H2.2I give a short review of notations and basics of polynomials and Grobner bases; 
Section |2~31 reviews the basic ideas of F5. 

For a more detailed introduction on non-F5 basics we refer the reader to [19j . Readers familiar 
with these topics may want to skim this section for notation and terminology. 

2.1. Polynomial basics. Let K. be a field, V :— K\x\ the polynomial ring over K. in the variables 
x := (xi, . . . , x n ). Let T denote the set of terms {x a } C V, where x a := Y\i=i x ? i an< ^ Q ' ; ^- 

A polynomial p over /C is a finite /C-linear combination of terms, i.e. p — ^ Q a a x a eP,o n G K,. 
The degree of p is the integer deg(p) = max{«i + ■ • ■ + a n a a =/= 0} for p =/= and deg(p) = —1 for 
p = 0. 

In this paper > denotes a fixed admissible ordering on the terms T. W.r.t. > we can write any 
nonzero p in a unique way as 

p = a a x a + apx 13 + . . . + a 1 x 1 , x a > x 13 > ■ ■ ■ > x 1 

where a a ,ap, . . . ,a 7 G /C\{0}. We define the head term of p HT(p) = x a and the head coefficient 
of p UC{p) = a a . 

2.2. Grobner basics. We work with homogeneous ideals / in V . For any S C V let HT(S) := 
(HT(p) | p G S\{0}}. A finite set G is called a Grobner basis of an ideal I if G C / and HT(I) = 
HT(G). Let p G V. If p = or there exist A, G V,qi G G such that p — 5^ =1 A*<ft and 
HT(p) > HT(A i q.j) for all nonzero qi, then we say that there exists a standard representation of p 
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w.r.t. G, or that p has a standard representation w.r.t G. We generally omit the phrase "w.r.t. G" 
when it is clear from the context. 

Let Pi,Pj G V. We define the s-polynomial of the critical pair (pi,pj) to be 

where 7tf := 1cm (HT(pi), HTfo)). 



Theorem 1. Let I be an ideal in V and G C I finite. G is a Grobner basis of I iff for allpi,pj € G 
Pij has a standard representation. 

Proof. See Theorem 5.64 and Corollary 5.65 in [1 pp. 219-221]. □ 

In addition to inventing the first algorithm to compute Grobner bases, Buchberger discovered 
two relatively efficient criteria that imply when one can skip an s-polynomial reduction [7,9|. We 
will refer occasionally to the second of these criteria. 

Theorem 2 (Buchberger's 1cm criterion). Let G C V be finite, and Pi,Pj,Pk € P. If 

(A) HT(p fe ) | fcm(HT( K ),HT( Pj )), and 

(B) pik and pjk have standard representations w.r.t. G, 
then p^ also has a standard representation w.r.t. G. 

In the homogeneous case one can define a d- Grobner basis Gd of an ideal /: This is a basis of 
I for which all s-polynomials up to degree d have standard representations (cf. Definition 10.40 
in [1 p. 473]). 

The following definition is crucial for understanding the problem of termination of F5. 

Definition 3. Let G be a finite set of polynomials in V . We say that p € G is redundant if there 
exists an element p' G G such that p' ^ p and HT(p') | HT(p). 

Remark 4. While computing a Grobner basis, a Buchberger-style algorithm does not add poly- 
nomials that are redundant at the moment they are added to the basis, although the addition of 
other polynomials to the basis later on may render them redundant. This ensures termination, as 
it expands the ideal of leading monomials, and V is Noetherian. However, F5 adds many elements 
that are redundant even when they are added to the basis; see Section [3l 

It is easy and effective to interreduce the elements of the initial ideal before F5 starts, so that 
the input contains only non-redundant polynomials; in all that follows, we assume that this is the 
case. However, even this does not prevent F5 from generating redundant polynomials. 

Finally, we denote by ip(p, G) the normal form of p with respect to the Grobner basis G. 

2.3. F5 basics. It is beyond the scope of this paper to delve into all the details of F5; for a more 
detailed discussion we refer the reader to [IB], [H], and [T3]. In particular, we do not consider the 
details of correctness for F5, which are addressed from two different perspectives in |16| and |13| . 
This paper is concerned with showing that the algorithm can be modified so that termination is 
guaranteed, and that the modification does not disrupt the correctness of the algorithm. 

In order to make the explanations more focused and concise, we now adapt some basic definitions 
and notation of [TBJ. Let Fj be the i-th canonical generator of V m . Denote T = U^Tj where 

= {tFi | t € T} and R = T x V. Define <, the extension of < to T, by tFi -< uFj iff 
(1) i > j, or 
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(2) i = j and t < u. 

It is easy to show that -< is a well-ordering of T, which implies that there exists a minimal repre- 
sentation in terms of the generators. 

Definition 5. Let p £ V and t £ T. We say that tFi is the signature of p if there exist hi, ... , h m £ 
V such that each of the following holds: 

• V = J2k=i hkfk and HT(fti) = t, and 

• for any H jt . . . ,H m £V such that p = YJk=j H kfk and Hj ^ 0, we have tF { < WT(Hj)Fj. 
Definition 6. Borrowing from [24], we call the element 

r = (tFi,p) £ R 

of |16] a labeled polynomial. (It is referred to as the representation of a polynomial in [16J.) We also 
denote 

(1) the polynomial part of r poly(r) = p, 

(2) the signature of r S(r) = tFi, and 

(3) the signature term of r ST(r) = t, and 

(4) the index of r index(r) = i. 

Following [TSj , we extend the following operators to R: 

(1) HT(r) = HT(p). 

(2) HC(r) = HC(p). 

(3) deg(r) = deg(p). 

Let / c £ K, \ £ T, r = (tFi,p) £ R. Then we define the following operations on R resp. T: 

(1) cr={tFi,cp), 

(2) Xr = (XtFi, Xp), 

(3) X(tF t ) = (At)Fj. 

Caveat lector: Although we call (S(r) the signature of r in Definition El it might not be the 
signature of poly(r) as defined in Definition[SJ If the input is non-regular, it can happen (and does) 
that F5 reduces an s-polynomial to zero. The reductions are all with respect to lower signatures, 
so we have 

#G 

polyfaj) = ^2 h k po\y(r k ) 
fc=i 

where h k ^ implies that S(WT(hk) ■ r k ) -< S(rij). The signature of nj is thus no larger than 
maxh k jto{S(WT(hk) ■ fk)}\ that is, the signature of is strictly smaller than S(rij). 

On the other hand, Propositions 171 and ITUl show that the algorithm does try to ensure that S(r) 
is the signature of poly(r). The proof of Proposition [7] is evident from inspection of the algorithm. 

Proposition 7. Let the list F = (/i, . . . , f m ) £ V m be the input of F5. For any labeled polynomial 
r = (tFi,p), t £ T, 1 < i < m, computed by the algorithm, there exist hi, ... , h m £ V such that 

(1) p = hxfi + ■ • • + h rn f m , 

(2) h\ = . . . = hi-i = 0, and 

(3) ST(r) = HT(^) = t. 

Let G — {ti , . . . , r na } C V. We denote poly(G) = {poly(ri), . . . , poly(r„ G )}. 
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Definition 8. Let r, r\, . . . , r„ G € i?, G — {r\ , . . . , r„ G }. Assume poly(r) 7^ 0. We say that r has 
a standard representation w.r.t. G if there exist Ai, . . . , A„ G <E V such that 

n G 



HT(r) > HT(A l )HT(r l ) for all i, and 5(r) >- HT(A i )5(r i ) for all z except possibly one, say i , 
where S(r) — S(ri ) and Xt = 1. We generally omit the phrase "w.r.t. G" when it is clear from 
the context. 

Remark 9. The standard representation of a labeled polynomial r has two properties: 

(1) The polynomial part of r has a standard representation as defined in Section 12.21 and 

(2) the signatures of the multiples of the are not greater than the signature of r. 

This second property makes the standard representation of a labeled polynomial more restrictive 
than that of a polynomial. 

Proposition 10. Let the list F = (/1, . . . , f m ) € V m be the input of F5. For any labeled polynomial 
r that is computed by the algorithm, if r does not have a standard representation w.r.t. G, then S(r) 
is the signature of poly (r). 

In other words, even if S(r) is not the signature of poly(r), the only time this can happen is 
when r already has a standard representation, so it need not be computed. On the other hand, the 
converse is false: once the algorithm ceases to reduce poly(r), r does have a standard representation, 
and S(r) remains the signature of poly(r). 

Proof. We show the contrapositive. Suppose that there exists some r £ G such that S(r) is not the 
signature of p = poly(r). Of all the r satisfying this property, choose one such that S(r) is minimal. 
Suppose S(r) = tFi. By hypothesis, we can find hj, . . . , h m € V such that 



Is hkfk a standard representation of r w.r.t. Gl Probably not, but it is clear that for any k 
such that HT(fr fc )HT(/fc) > HT(p), there exists £ such that HT(/i fc )HT(/ fc ) = HT(^)HT(/ £ ). The 
signature of the corresponding s-polynomial pki is obviously smaller than tFi, so the hypothesis 
that 5(r) is minimal, along with inspection of the algorithm, implies that r^i has a standard 
representation w.r.t. G. Proceeding in this manner, we can rewrite ^ hkfk repeatedly until we 
have a standard representation of r w.r.t. G. □ 

Definition 11. Let r l = (tiF k) pi),rj = (tjFi,pj) E R. If ^\ j ri) tiF k ^ H ^' ( J r , tjF t , then we define 
the s-polynomial of ri and rj by := (m! ,Pij) where 



and 7ij - = lcm(HT(r 4 ),HT(r J )). 

All polynomials are kept monic in F5; thus we always assume in the following that HC(pi) = 
HC(pj) = 1 for pi ^ 9^ Pj. Moreover we always assume 7y to denote the least common multiple 
of the head terms of the two considered polynomial parts used to compute rij . 

Next we review the two criteria used in F5 to reject critical pairs which are not needed for further 
computations. 




1 n 




hj ^ 0, and i < j or [i — j and t > WT(hj)] . 



k =J 
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Definition 12. Let G — {ri, . . . ,r nG } be a set of labeled polynomials, and Uk £ T. We say that 
MfcTfc is detected by Faugere's Criterion if there exists r £ G such that 

(1) index(r) > index(r^) and 

(2) HT(r) | u k ST(r k ). 

Definition 13. Let G — {ri, . . . ,r na } be a set of labeled polynomials, and £ T. We say that 
MfcTfc is detected by the Rewritten Criterion if there exists r a £ G such that 

(1) mdex(r a ) = index(r fc ), 

(2) a > k, and 

(3) ST(r Q ) | u k ST(r k ). 

Next we can give the main theorem for the idea of F5. Recall that we consider only homogeneous 
ideals. 

Theorem 14. Let I = {fi, ■ ■ . , f m ) be an ideal in V, and G = {r l7 ...,r„ G } a set of labeled 
polynomials generated by the F5 algorithm (in that order) such that fi £ poly(G) for 1 < i < m. 
Let d £ N. Suppose that for any pair r i: rj such that degr.jj < d and r^ — u^i — UjTj, one of the 
following holds: 

(1) Ufcrfc is detected by Faugere's Criterion for some k £ {i,j}, 

(2) Mfcrfc is detected by the Rewritten Criterion for some k £ {i,j}, or 

(3) rij has a standard representation. 
Then poly(G) is a d-Grdbner basis of I. 

Proof. See Theorem 1 in [16], Theorem 3.4.2 in [T7] and Theorem 21 in [13]. □ 
Remark 15. 

(1) Requiring a standard representation of a labeled polynomial is stricter than the criterion of 
Theorem [TJ but when used carefully, any computational penalty imposed by this stronger 
condition is negligible when compared to the benefit from the two criteria it enables. 

(2) It is possible that does not have a standard representation (cf. Proposition 17 in |13| ) 
at the time either Criterion rejects (r,,rj). Since F5 computes the elements degree- by- 
degree, computations of the current degree add new elements such that r.^ has a standard 
representation w.r.t. the current Grobner basis poly(G) before the next degree step is 
computed. Thus, at the end of each such step, we have computed a d-Grobner basis of I. 

Next we give a small example which shows how the criteria work during the computation of a 
Grobner basis in F5. 

Example 16. Let > be the degree reverse lexicographical ordering with x > y > z on Q[x,y,z]. 
Let / be the ideal generated by the following three polynomials: 

pi = xyz - y 2 z, 

P2 = x 2 - yz, 

P3 = V 2 ~ xz. 

Let the corresponding labeled polynomials be r, = (Fj,j>j). For the input F — (pi,P2,P3), F5 
computes a Grobner basis of (p2,P3) a s a first step: Since ST^^) = y 2 = HT(rs), ^3 is discarded 
by Faugere's Criterion. Thus {^2,^3} is already a Grobner basis of (^2,^3)- 

Next the Grobner basis of / is computed, i.e. r\ enters the algorithm: Computing ^,3 we get 
a new element: r^ = (yFi,xz 3 — yz 3 ). r± 2 is n °t discarded by any criterion, but reduces to zero. 
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Nevertheless its signature is recorded^ thus we still have S(ri^) = xF\ stored in the list of rules 
to check subsequent elements. 

Next check all s-polynomials with r± sorted by increasing signature: 

(1) Since S{r^ i) = y 2 F 1 , r^i is discarded by Faugere's Criterion using HT(r 3 ) = y 2 . 

(2) Since 5(r4 2 ) = xyFi, r± 2 is discarded by the Rewritten Criterion due to 5(r 12 ) = xF-l, 
r\ 2 being computed after r^. 

(3) Since £(74 3) = y 3 Fi, r4,3 is discarded by Faugere's Criterion using HT(r3) = y 2 . 

The algorithm now concludes with G — {7*1, T2, T3, r^] where poly(G) is a Grobner basis of I. 

3. Analysis of the problem 

The root of the problem lies in the algorithm's reduction subalgorithms, so Section |3~T1 reviews 
these in detail. In Section l3~2| we show how the criteria force the reduction algorithms not only 
to add redundant polynomials to the basis, but to do so in a way that does not expand the ideal 
of leading monomials (Example I18[) ! One might try to modify the algorithm by simply discarding 
redundant polynomials, but Section 13.31 shows that this breaks the algorithm's correctness. This 
analysis will subsequently provide insights on how to solve the problem. 

Throughout this section, let the set of labeled polynomials computed by F5 at a given moment 
be denoted G — {r± , . . . , r„ G } . 

3.1. F5's reduction algorithm. For convenience, let us summarize the reduction subalgorithms 
in some detail here. Let i be the current iteration index of F5. All newly computed labeled 
polynomials r satisfy indcx(r) = i. Let Gi+\ denote the set of elements of G with index > i. We 
are interested in Reduction, TopReduction and IsReducible. F5 sorts s-polynomials by degree, 
and supplies to Reduction a set F of s-polynomials of minimal degree d. Let r G F. 

(1) First, Reduction replaces the polynomial part of r with its normal form with respect to 
Gi+i. This clearly does not affect the property S(r) — ST(r)Fj. Reduction then invokes 
TopReduction on r. 

(2) TopReduction reduces poly(r) w.r.t. Gi, but invokes IsReducible to identify reducers. 
TopReduction terminates whenever poly(r) = or IsReducible finds no suitable reducers. 

(3) IsReducible checks all elements r rc d S G such that index(r roc i) = i. 

(a) If there exists u re( j £ T such that u rcc iHT(r rcc i) = HT(r) then ii re d<S(r re d) is checked 
by both Faugere's Criterion and the Rewritten Criterion. 

(a) If neither criterion holds, the reduction takes place, but a further check is neces- 
sary to preserve S{r) = ST(r)IV If S{r) >~ u re d<S(r re d)) then it rewrites poly(r): 

r = (,S(r),poly(r) - u re dPoly(r red )). 

If S(r) -< u Ie dS(r Te d), then r is not changed, but a new labeled polynomial is 
computed and added to F for further reductions, 

r' = (urcdS(ncd),w rc dPoly(r red ) -poly(r)). 

The algorithm adds S(r') to the list of rules and continues with r. 

If Ured^red is detected by one of the criteria, then the reduction does not take 

place, and the search for a reducer continues. 

(b) If there is no possible reducer left to be checked then r is added to G if poly(r) 7^ 0. 



Failing to record the signature of a polynomial reduced to zero is an implementation error that can lead to an 
infinite loop. 
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Note that if S(r) = u re dS(r le d) then w re d^red is rewritable by r, thus Case (3)(a)(/3) avoids this 
situation. 

3.2. What is the problem with termination? The difficulty with termination arises from Case 
(3) (a) 08) above. 

Situation 17. Recall that Rd is the set of labeled polynomials returned by Reduction and added 
to G. Suppose that Rd 7^ and for every element r £ Rd, HT(poly(r)) is in the ideal generated by 
HT(poly(G)). 

Example 18. Situation[T7]is not a mere hypothetical: as described in Section 3.5 of [17], an example 
appears in Section 8 of |16| . which computes a Grobner basis of (yz 3 — x 2 t 2 ,xz 2 — y 2 t,x 2 y — z 2 t). 
Without repeating the details, at degree 7, F5 adds r§ to G, with HT(rg) = y 5 t 2 . At degree 8, 
however, Reduction returns R$ = {rio}, with HT(rio) = y 6 t 2 . This is due to the fact that the 
reduction of rio by yr% is rejected by the algorithm's criteria, and the reduction does not take place. 
In other words, rio is added to G even though poly(rio) is redundant in poly(G). 

Definition 19. A labeled polynomial r computed in F5 is called redundant if, when Reduction 
returns r, we have poly(r) redundant w.r.t. poly(G). 

Lemma 20. // Rd satisfies Situation \17\ and r £ Rd, then we can find r^ £ G such that r% is not 
redundant in G and HT(rt) | HT(r). 

Proof. If a reducer rj of r is redundant, then there has to exist another element r^ such that 
HT(rfc) | HT(rj) and thus HT(rfc) | HT(r). Follow this chain of divisibility down to the minimal 
degree; we need to show that there do not exist two polynomials rj. r^ of minimal degree such 
that HT(rj) = HT(rfc). Assume to the contrary that there exist rj,rk £ G of minimal degree such 
that HT(rj) = HT(rfc). Clearly, the reduction of one by the other in IsReducible was forbidden; 
without loss of generality, we may assume that was computed before rj , so the reduction of rj 
by r,t was forbidden. There are three possibilities: 

(1) If index(rj.) > index(rj), to the contrary, IsReducible cannot interfere with this reduc- 
tion, because such reductions are always carried out by the normal form computation in 
Reduction. 

(2) If S(rk) is rejected by the Rewritten Criterion, then there exists r' such that ST(r') | ST(rfc), 
and r' was computed after rfe. (That r' was computed after r^ follows from Definition I13[ 
where a > k.) As F5 computes incrementally on the degree and ST(r') | ST(r-fc), it follows 
that deg(r') = deg(rfc). Hence ST(r') = ST(rfc). Thus the Rewritten Criterion would have 
rejected the computation of r' , again a contradiction. 

(3) If iS(rfc) is rejected by Faugere's Criterion, to the contrary, r^ should not have been com- 
puted in the first place. 

Thus HT(rj) ^ HT(rt). It follows that we arrive at a non-redundant reducer after finitely many 
steps. □ 

Lemma 21. Denote by Rd the result of Reduction at degree d. There exists m £ N and an input 
F = (fx, . . . , f m ) and a degree d such that if poly(G) is a (d — 1)-Grdbner basis of (fi, . . . , f m ), 
then 

(A) R d ^ 0, and 

(B) i?T(poly(G U R d )) = #T(poly(G)). 
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Proof. Such an input F is given in Example 1181 once reduction concludes for d = 8, HT(rs) | 
HT(no), so HT(poly(G)) = HT(poly(G U R 8 )). □ 

Remark 22. In [TBI Corollary 2], it is argued that termination of F5 follows from the (unproved) 
assertion that for any d, if no polynomial is reduced to zero, then HT(poly(G)) ^ HT(poly(GUi?d)). 
But in Example 1181 HT(poly(G)) = HT(poly(GUi?8)); even though there was no reduction to zero! 
Thus, Theorem 2 (and, by extension, Corollary 2) of [TFJ is incorrect: termination of F5 is unproved, 
even for regular sequences, as there could be infinitely many steps where new redundant polynomials 
are added to G. By contrast, a Buchberger-style algorithm always expands the monomial ideal when 
a polynomial does not reduce to zero; this ensures its termination. 

Having shown that there is a problem with the proof of termination, we can now turn our 
attention to devising a solution. 

3.3. To sort the wheat from the chaff . . . isn't that easy! The failure of F5 to expand the 
ideal of leading monomials raises the possibility of an infinite loop of redundant labeled polynomials. 
However, we cannot ignore them. 

Example 23. Suppose we modify the algorithm to discard critical pairs with at least one redundant 
labeled polynomial. Consider a polynomial ring in a field of characteristic 7583. 

(1) For Katsura-5, the algorithm no longer terminates, but computes an increasing list of 
polynomials with head terms x\x± with signatures xix\x$x§ for k > 1. 

(2) For Cyclic-8, the algorithm terminates, but its output is not a Grobner basis! 

How can critical pairs involving "redundant" polynomials can be necessary? 

Definition 24. A critical pair (r^rj) is a GB-critical pair if neither rj nor rj is redundant. If a 
critical pair is not a GB-critical pair, then we call it an F5-critical pair. 

We now come to the main theoretical result of this paper. 

Theorem 25. If (r,,rj) is an F5-critical pair, then one of the following statements holds at the 
moment of creation of rij : 

(A) poly (rij) already has a standard representation. 

(B) There exists a GB-critical pair (ru,ri), a set W C {1, . . . ,uq\ , and terms X w (for all 
w G W) such that 

poly(ry) = poly(r w ) + ^ X w poly(r w ), (3.1) 

w 

7ij = 7k£ and j ke > \ w WT(r w ) for all w. 

Theorem [25] implies that an F5-critical pair might not generate a redundant polynomial: it might 
rewrite a GB-critical pair which is not computed. Suppose, for example, that the algorithm adds rj 
to G, where r% is redundant with r k G G, perhaps because for u G T such that uHT(r^) = HT(rj), 
we have S(u • r k ) >- <S(?"i). In this case, F5 will generate a new, reduced polynomial with the 
larger signature; since the new polynomial has signature S(u ■ r^), the Rewritten Criterion will 
subsequently reject u ■ r^. It is not uncommon that the algorithm later encounters some rg G G 
where r^ is necessary for the Grobner basis, but HT(r^) divides j^e- In this case, the Rewritten 
Criterion forbids the algorithm from computing r k e, yet we can compute ru. In terms of the 
Macaulay matrix [T6l[20l[2T] . the algorithm selects the row corresponding to H ^ ; r j instead of 
the row corresponding to Due to this choice, the notions of "redundant" and "necessary" 
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critical pairs are somewhat ambiguous in F5: is fj necessary to satisfy the properties of a Grobner 
basis, or to ensure correctness of the algorithm? On the other hand, the notions of F5- and GB- 
critical pairs are absolute. 

To prove Theorem I25[ we need the following observation: 

Lemma 26. Let r^rj £ G computed by F5, and assume that WT(rj) | HT(r,). Then Spol does 
not generate an s-polynomial for fa,?^). 

Proof. We have assumed that the input is interreduced, so poly(ri) is not in the input. Since 
HTfa) I HT(r-j) there exists u £ T such that uKT(r ) = HT(r;). Since the reduction of poly(ri) 
by upoly(rj) was rejected, uS(rj) was detected by one of the criteria. It will be detected again in 
CritPair or Spol. Thus Spol will not generate ry. □ 

Proof of Theorem \25l Assume that fj and Vj are both redundant; the case where only (resp. rj) 
is redundant is similar. By Lemma 1201 there exists for rj (resp. rj) at least one non-redundant 
reducer r k (resp. rg). By Lemma \26\ we may assume that r, and rj are of degree smaller than ry. 
Using the fact that poly(G) is a e?-Grobner basis for d = max(degrj,degrj), we can write 

poly(r;) = Aifcpoly(r fe ) + ^ A„poly(r u ) 

ueu 

P°ly fa) = VP ol yfa) + ^2 X v po\y(r v ), 

such that 

HT(ri) = \ lk UT(r k ) > A„HT(r„) and 
HTfa) = \ je H.T{r e ) > A„HT(r„) 

where U, V C {1, . . . , uq\. As 7^ | 7^, the representations of poly(rj) and poly(rj) above imply 
that there exists A £ T such that 

p 01 ^ = HT^ poly(ri) " fflfe poly( ^ 

= Apoly (r k e) + ^ A w poly(r„,) (3.2) 

where W = U U V and X w — jj^ry A u for w £ U\V. X w = HT ^ J r ^ A„ for 10 £ V\f7, and A M = 
H t'^.) A m — jrr^rjA^ for w £ C/ n V. In Equation (|3.2[) we have to distinguish two cases: 

(1) If A > 1 then deg(r^) < deg(ry), thus r k e is already computed (or rewritten) using a lower 
degree computation, which has already finished. It follows that there exists a standard 
representation of poly(r^) and thus a standard representation of poly(r.y ■). 

(2) If A = 1 then (A) holds if poly(rfci) is already computed by F5; otherwise (B) holds. 

□ 

We can now explain why discarding redundant polynomials wreaks havoc in the algorithm. 

Situation 27. Let fa,fj) be an F5-critical pair. Suppose that all GB-critical pairs (r k ,ri) cor- 
responding to case (B) of Theorem \25\ are rejected by one of F5's criteria, but lack a standard 
representation. 

Situation 1271 is possible if, for example, the Rewritten Criterion rejects all the (r k ,r^). 
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Corollary 28. In Situation \27\ it is necessary for the correctness of F5 to compute a standard 
representation ofr^. 

Proof. Since po\y(rke) lacks a standard representation, and the algorithm's criteria have rejected the 
pair (rk,rg), then it is necessary to compute a standard representation of r^. Once the algorithm 
does so, we can rewrite (|3.1[) to obtain a standard representation of poly(rj^). □ 

In other words, "redundant" polynomials are necessary in F5. 

4. Variants that ensure termination 

Since we cannot rely on an expanding monomial ideal, a different approach to ensure termination 
could be to set or compute a degree bound. Since a Grobner basis is finite, its elements have a 
maximal degree. Correspondingly, there exists a maximal possible degree cZgb of a critical pair that 
generates a necessary polynomial. Once we complete degree c?qb, no new, non- redundant data for 
the Grobner basis would be computed from the remaining pairs, so we can terminate the algorithm. 
The problem lies with identifying g?gb, which is rarely known beforehand, if ever0 

Before describing the new variant that follows from these ideas above, we should review two 
known approaches, along with some drawbacks of each. 

4.1. F5t: Reduction to zero. In [17J, Gash suggests the following approach, which re-introduces 
a limited amount of reduction to zero. Once the degree of the polynomials exceeds 2M, where M 
is the Macaulay bound for regular sequences [H[M|) start storing redundant polynomials in a set 
D. Whenever subalgorithm Reduction returns a nonempty set Rd that does not expand the ideal 
of leading monomials, reduce all elements of Rd completely w.r.t. G U D and store any non-zero 
results in D instead of adding them to G. Since complete reduction can destroy the relationship 
between a polynomial and its signature, the rewrite rules that correspond to them are also deleted. 
Subsequently, s-polynomials built using an element of D are reduced without regard to criterion, 
and those that do not reduce to zero are also added to D, generating new critical pairs. Gash called 
the resulting variant F5t. 

One can identify four drawbacks of this approach: 

(1) The re-introduction of zero-reductions incurs a performance penalty. In Gash's experiments, 
this penalty was minimal, but these were performed on relatively small systems without 
many redundant polynomials. In some systems, such as Katsura-9, F5 works with hundreds 
of redundant polynomials. 

(2) It keeps track of two different lists for generating critical pairs and uses a completely new 
reduction process. An implementation must add a significant amount of complicated code 
beyond the original F5 algorithm. 

(3) It has to abandon some signatures due to the new, signature-corrupting reduction process. 
Thus, a large number of unnecessary critical pairs can be considered. 

(4) The use of 2M to control the size of D is an imprecise, ad-hoc patch. In some experiments 
from p.7j, F5t terminated on its own before polynomials reached degree 2M; for other input 
systems, F5t yielded polynomials well beyond the 2M bound, and a higher bound would 
have been desirable. 



2 Another algorithm that computes a degree bound is MXL3 1221 , but its mechanism is designed for zero-dimensional 
systems over a field of characteristic 2. It is not appropriate for the general case, whereas the approaches that we 
study here are. 
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4.2. F5B: Use Buchberger's 1cm criterion. In pQ, Ars suggests using Buchberger's 1cm criterion 
to determine a degree bound. 

• Initialize a global variable ds — storing a degree. 

• Keep a second list of critical pairs, P*, used only to determine a degree bound. 

• When adding new elements to G, store a copy of each critical pair not detected by Buch- 
berger's 1cm criterion in P* . Remove any previously-stored pairs that are detected by 
Buchberger's 1cm criterion, and store the highest degree of an element of P* in ds- 

If the degree of all critical pairs in P exceeds ds, then a straightforward application of Buchberger's 
lcm criterion implies that the algorithm has computed a Grobner basis, so it can terminate. We 
call this variant F5B. 

It is important to maintain the distinction between the two lists of critical pairs. Otherwise, the 
correctness of the algorithm is no longer assured: Buchberger's criteria ignore the signatures, so P* 
lacks elements needed on account of Situation [27] 

While elegant, this approach has one clear drawback. Every critical pair is computed and checked 
twice: once for Buchberger's lcm criterion, and again for the F5 criteria. Although Faugere's Cri- 
terion also checks for divisibility, it checks only polynomials of smaller index, whereas Buchberger's 
criterion checks all polynomials, and in most systems the number of polynomials of equal index is 
much larger than the total of all polynomials having lower index. Indeed, we will see in Section 
14.51 that this seemingly innocuous check can accumulate a significant time penalty. This would be 
acceptable if the algorithm routinely used ds to terminate, but F5 generally terminates from its 
own internal mechanisms before d = ds! Thus, except for pathological cases, the penalty for this 
short-circuiting mechanism is not compensated by a discernible benefit. 

4.3. F5+: Use F5's criteria on non-redundant critical pairs. We now describe a variant 
that uses information from F5 itself, along with the theory developed in Section [3j to reduce, if not 
eliminate, the penalty necessary to force termination. We restate only those algorithms of [16J that 
differ from the original (and the differences are in fact minor). 

The fundamental motivation of this approach stems from the fact that a polynomial is redundant 
if and only if TopReduction rejects a reductor on account of one of the F5 criteria. Understood 
correctly, this means that F5 "knows" at this point whether a polynomial is redundant. We would 
like to ensure that it does not "forget" this fact. As long as this information remains available to 
the algorithm, identifying GB- and F5-critical pairs will be trivial. Thus, our tasks are: 

(1) Modify the data structures to flag a labeled polynomial as redundant or non-redundant. 

(2) Use this flag to distinguish F5- and GB-critical pairs. 

(3) Use the GB-critical pairs to decide when to terminate. 

We address each of these in turn. 

To distinguish between redundant and non-redundant labeled polynomials, we add a third, 
boolean field to the structure of a labeled polynomial. We mark a redundant labeled polyno- 
mial with 6 = 1, and a non-redundant one with 6 = 0. Without loss of generality, the inputs are 
non-redundant, so the first line of subalgorithm F5 can change to 

n := (Fi./i.O) Gi?x {0,1} 
For all other labeled polynomials, the value of 6 is set to in algorithm Spol, then defined by the 
behaviour of the Reduction subalgorithm; see below. 

The next step is to detect redundant polynomials; we do this in IsReducible. In an unmodified 
F5, the return value of IsReducible is either a labeled polynomial fj. (a polynomial that reduces 
r) or 0. The return value can have two meanings: 



MODIFYING FAUGERE'S F5 ALGORITHM TO ENSURE TERMINATION 



13 



(1) There exists no reducer of the input. 

(2) There exist reducers of the input, but their reductions are rejected. 

Algorithm [TJ which replaces the original IsReducible subalgorithm, distinguishes these two possi- 
bilities by adding a boolean to the output: b — in case ([T]) and 6=1 otherwise. We also need to 
modify subalgorithm TopReduction to use this new data; see Algorithm [2] 

We now describe the main routine of the new variant, which fulfills the following conditions: 

(1) Compute as low a degree bound as possible. 

(2) Minimize any penalty to the algorithm's performance. 

An easy way to estimate do would be to compute the highest degree of a GB-critical pair. Although 
this would be correct, experience suggests that, in general, it is much higher than necessary (see 
Table [1] in Section |4~51) . Instead, the new variant will use the criteria of the F5 algorithm to identify 
GB-critical pairs that probably reduce to zero. How can we identify such pairs? The following 
method seems intuitively correct: when all GB-critical pairs are rejected by one of the F5 criteria. 

However, Situation 1271 implies that this intuition may be incorrect. Thus, once the algorithm 
reaches that degree (and not earlier), it uses Buchberger's 1cm criterion to decide whether the 
remaining GB-critical pairs reduce to zero. If it can verify this, then the algorithm can terminate. 

This differs from the approach of pQ in two important ways. 

(1) Rather than checking all pairs against the 1cm criterion, it checks only GB-critical pairs 
that F5 also rejects as unnecessary. After all, it follows from Theorem [231 that F5-critical 
pairs can be necessary only if they substitute for a GB-critical pair. 

(2) It checks the GB-critical pairs only once the F5 criteria suggest that it should terminate. 

We call this variant F5+; see Algorithm [3l 
Algorithm 1 IsReducible 



b := 1 
return (0, b) 



Remark 29. An implementation of F5+ has to take care when checking Buchberger's 1cm criterion, 
on account of the phenomenon of Buchberger triples [31 p. 229]. In [T], this is implemented similarly 
to the "Update" algorithm of [51ll8|. The current F5+ takes a more traditional route; it records all 
critical pairs that have generated s-polynomials. The burden on memory is minimal. 

4.4. Correctness and termination of F5+. As a last step we have to show that F5+ terminates 
correctly. 




ri , a labeled polynomial of R 



b:=0 



for j from 1 to r do 



if (u := g T ) then 

if (neither criterion detects (n ,n )) then 

return (r; . , 0) 
else 
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Algorithm 2 TopReduction 



Input 



Tk , a labeled polynomial of R 
G, a list of elements of R 

ken 

if. a normal form 
if poly(r/c ) = then 

return (0,0) 
(r',b) :=IsReducible(rfc , G, k, if) 
if r' = then 

r ko := (S (r ko ) , ^ poly (r fco ) , b) 

return (r feo ,0) 
else 

rki = r' 

HT(r t „) 
"" HT(r fcl ) 

if uS(r kl ) -< <S(r feo ) then 

r kg := (S(r fe J,poly(r feo ) -upoly(r fel ),6) 

return (0, {r fco }) 
else 

TV := N + 1 

r N := (u,S(r fcl ),Mpoly(r fcl ) - poly(r feo ), b) 
Add Rule (rjv) 
return (0, {r N , ru }) 



Theorem 30. If F5+ terminates, the result is a Grobner basis of the input. 

Proof. This follows from Buchberger's 1cm criterion. □ 

Theorem 31. For a given homogeneous ideal I as input, F5+ terminates after finitely many steps. 

Proof. We first claim that after generating new critical pairs for P in lines [22H291 F5+ satisfies 
ffP < oo at line [30l and thus satisfies #P < oo when the loop at line [14] iterates anew. To show 
this, we will show that at any given degree d, the algorithm generates only finitely many polynomials 
and critical pairs. We proceed by induction on d; certainly #P < oo after the loop in lines [THT^I 
Assume therefore that #P < oo at linefMl By the assumption that #P < oo, we have #Pd < oo, 
so Spol generates only finitely many new polynomials. We now consider Reduction; let r 6 ToDo. 

(1) If poly(r) = 0, then r is effectively discarded; the algorithm does not add it to G, nor use 
it to generate new critical pairs. 

(2) If poly(r) 7^ 0, then IsReducible checks for possible reducers: 

(a) If no reducer is returned, then r is returned and added to G. All newly computed 
critical pairs generated by r have degree > d; their number is finite because G is 
currently finite. 

(b) If IsReducible returns r rc d e G such that there exists w ro d £ T satisfying w ro dHT(r rcc i) = 
HT(r) and S(r) >- u Te dS(r re d), then poly(r) — u r cdPoly(r rc d) replaces poly(r) in r, and 

r is checked for further reductions. Note that HT(poly(r)) has decreased. 
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Algorithm 3 F5+ 



i g N 

Input: ^ fi G K,[x\ 

Gi+i C R x /C[x], such that poly(Gj + i) is a Grobner basis of Id(/j + i, . . . , f m ) 
n := (Fj./j.O) 

:=NF(.,poly(G i+1 )) 
G, := G i+1 U {rj 

{P is the usual set of pairs; P* is the set of GB-pairs detected by the F5 criterion} 
P := 
P* := 

for rj G G, + i do 

p := CritPair(rj, rj, i, <Pi+i) 
if p = and rj non-redundant then 

Add (lcm(HT(poly(r J )),HT(poly(r J ))),r 4 ,r J ) to P* 
else 

Add p to P 
Sort P by degree 
while P ^ do 
d := deg(first(P)) 

Discard from P* all pairs that are not of maximal degree 

if d < max{deg(p) : p G P*} or Bp G P* that does not satisfy Buchberger's 1cm criterion 
then 

P d :={peP: deg(p) - d} 

P := P\Pd 
P := Spol(Pi) 

Pd := Reduction(P, Gj, i, ^,+i) 
for r G P<j do 
for rj G Gi do 
p := CritPair(r, rj,i, <fi+i) 
if p = and r, rj both non-redundant then 

Add (lcm(HT(poly(r)),HT(poly(rj))), r,r 3 ) to P* 
else 

Add p to P 
G, := G 4 U {r} 
Sort P by degree 
else 
P := 
return Gi 



(c) If IsReducible returns r ro d G G such that there exists w ro d G T satisfying u rc( jHT(r rc[ |) = 
HT(r) but Wrcd<S(r ro d) >~ S{r), then r is not changed, but kept for further reduc- 
tion checks. A new element r' — (u re dS(r le d), poly(r) — w re dpoly(r ro d)) is generated, 
and its signature <S(r') added to the list of rules. Note that degr' = degr and 
degST(r') =degST(r). 
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Only finitely many distinct reducers could lead to new elements r' . Since S(r') was 
added to the list of rules, the Rewritten Criterion implies that r ro( j will not be chosen 
again as a reducer of r. There are only finitely many signatures of degree d, so only 
finitely many new elements can be added in this way. 

It follows that in each degree step only finitely many new polynomials are computed, so only finitely 

many new critical pairs are generated. Hence #P < oo at line 1301 

To finish the proof we have to show that after finitely many steps, only F5-critical pairs are left 

in P. There can only be finitely many GB-critical pairs as their generating labeled polynomials 

have to be non-redundant. Since R is Noetherian, only finitely many non-redundant polynomials 

can be computed. 

Thus F5+ terminates after finitely many steps. □ 

4.5. Experimental results. We implemented these variants in the SINGULAR kernel to compare 
performance. (The F5 implementation in SINGULAR is still under development.) In Table Q] we 
compare timings and degree bounds for some examples. All systems are homogeneous and com- 
puted over a field of characteristic 32003. The random systems are generated using the function 
sparseHomogldeal from random, lib in SINGULAR; generating polynomials have a sparsity of 85- 
90% and degrees < 6. This data was recorded from a workstation running Gentoo Linux on an 
Intel® Xeon® X5460 CPU at 3.16GHz with 64 GB RAM. 

Table [1] shows that the tests for F5+ do not slow it down significantly. But this is expected, 
since the modifications add trivial overhead, and rely primarily on information that the algorithm 
already has available. 

The computed degrees in Table [1] bear some discussion. We have implemented F5+ in two 
different ways. Both are the same in that they estimate the maximum necessary degree by counting 
the maximal degree of a GB-critical pair not discarded by the CritPair subalgorithm. However, one 
can implement a slightly more efficient CritPair algorithm by discarding pairs that pass Faugere's 
Criterion, but are rewritable. (The basic F5 checks the Rewritten Criterion only in subalgorithm 
Spol.) Thus one might compute a different maximal degree of P* in each case: when CritPair 
discards only those pairs detected by Faugere's Criterion, we designate the maximal degree of P* 
as (If; when CritPair discards pairs detected by the Rewritten Criterion as well, we designate the 
maximal degree of P* as dpR. We denote the degree where the original F5 terminates by dps, and 
the maximal degree of a polynomial generated by rfmaxGB- Recall also that the maximal degree 
estimated by F5B is ds (Section 14. 2|) . 

It is always the case that rf max GB < ^fs; indeed, we will have <i m axGB < dA for any algorithm A 
that computes a Grobner basis of a homogeneous system incrementally by degree. 

On the other hand, it is always the case that dp, dFR < dp5] g^fs counts F5-critical pairs as well 
as GB-critical pairs, whereas dp,dpR count only GB-critical pairs that are not rejected by one or 
both of the F5 criteria. Thus F5+ always starts its manual check for termination no later than F5 
would terminate, and sometimes terminates before F5. For example, the termination mechanisms 
activate for F-855, Eco-10 and -11, and Cyclic-8, so F5B and F5+ both terminate at lower degree 
than F5. With little to no penalty, F5+ terminates first, but F5B terminates well after F5 in spite 
of the lower degree! Even in Katsura-n, where cf max GB = ds < dp = dFR = dp$, the termination 
mechanism of F5+ incurs almost no penalty, so its timings are equivalent to those of F5, whereas 
F5B is slower. In other examples, such as Cyclic-7 and (4,5,12), F5 and (therefore) F5+ terminate 
at or a little after the degree(s) predicted by dp and dFR, but before reaching the maximal degree 
computed by dg. 



Table 1. Timings (in seconds) & degrees of F5, F5B, and F5+ 



1 1 

h/xamples 


regular? 


F5 


F5B 


F5+ 


F5/F5B 


F5/F5+ 


j 2 

"maxGB 


d F 5 


, 4 

"GB-pair 


d B b 


d F b 


dFR 


Katsura-9 


yes 


39.95 


53.97 


40.23 


0.74 


0.99 


13 


16 


21 


13 


16 


16 


Katsura-10 


yes 


1,145.47 


1,407.92 


1,136.43 


0.80 


1.00 


15 


18 


26 


15 


18 


18 


F-855 


no 


9,831.81 


11,364.47 


9,793.17 


0.86 


1.00 


14 


18 


20 


17 


17 


16 


Eco-10 


no 


47.26 


57.97 


46.67 


0.82 


1.01 


15 


20 


23 


17 


17 


17 


Eco-11 


no 


1,117.13 


1,368.44 


1,072.47 


0.82 


1.04 


17 


23 


26 


19 


19 


19 


Cyclic-7 


no 


6.24 


9.18 


6.21 


0.67 


1.00 


19 


23 


28 


24 


23 


21 


Cyclic-8 


no 


3,791.54 


4,897.61 


3,772.66 


0.77 


1.00 


29 


34 


41 


33 


32 


30 


4,6,8 


no 


195.45 


204.88 


195.69 


0.95 


1.00 


22 


36 


42 


34 


34 


34 


5,4,8 


yes 


45.103 


46.930 


45.123 


0.96 


1.00 


20 


22 


35 


23 


20 


20 


6,4,8 


no 


46.180 


46.880 


46.247 


0.99 


1.00 


20 


20 


34 


22 


20 


20 


7,4,8 


no 


0.827 


0.780 


0.830 


1.060 


1.00 


14 


19 


27 


14 


17 


15 


8,3,8 


no 


122.972 


126.816 


123.000 


0.97 


1.00 


22 


37 


35 


26 


31 


29 


4,5,12 


no 


4.498 


5.680 


4.590 


0.79 


0.98 


29 


33 


37 


42 


32 


30 


6,5,12 


yes 


12.071 


21.150 


12.060 


0.57 


1.00 


50 


54 


73 


55 


54 


50 


8,4,12 


no 


46.122 


47.613 


47.750 


0.97 


0.97 


27 


35 


44 


30 


34 


29 


12,4,12 


no 


14.413 


14.897 


14.360 


0.97 


1.00 


42 


55 


60 


43 


53 


43 


4,3,16 


yes 


1.439 


1.403 


1.450 


1.03 


0.99 


15 


15 


23 


18 


15 


15 


6,3,16 


yes 


36.300 


37.136 


36.300 


0.98 


1.00 


10 


14 


23 


15 


14 


13 


8,3,16 


yes 


467.560 


471.737 


467.530 


0.99 


1.00 


12 


16 


21 


13 


15 


13 


12,3,16 


yes 


210.327 


206.441 


210.311 


1.02 


1.00 


21 


25 


34 


20 


24 


23 


4,3,20 


yes 


1.512 


1.680 


1.500 


0.90 


1.01 


16 


22 


24 


22 


21 


21 


6,4,20 


no 


1,142.433 


1,327.540 


1,144.370 


0.86 


1.00 


27 


37 


39 


29 


35 


31 


8,4,20 


no 


8.242 


8.230 


8.251 


1.00 


1.00 


35 


40 


48 


36 


40 


37 


12,3,20 


yes 


0.650 


0.693 


0.650 


0.94 


1.00 


22 


26 


34 


27 


26 


23 


16,3,20 


no 


2.054 


2.060 


2.050 


1.00 


1.00 


26 


26 


41 


27 


26 


26 



1 The notation (a, b, c) denotes a random system of a generators with maximal degree 6 in a polynomial ring of c variables. 

2 maximal degree in GB 

3 observed degree of termination of F5 

4 maximal degree of a GB-critical pair 

5 maximal degree estimated by Buchberger's 1cm criterion; see Section 14.21 

6 maximal degree of all GB-critical pairs not detected by Faugere's Criterion 

7 maximal degree of all GB-critical pairs not detected by Faugere's Criterion or the Rewritten Criterion 
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5. Concluding remarks, and a conjecture 

The new variant of F5 presented here is a straightforward solution to the problem of termination: 
it distinguishes F5- and GB-critical pairs and tracks the highest degree of a necessary GB-critical 
pair. Thus F5+ provides a self-generating, correct, and efficient termination mechanism in case F5 
does not terminate for some systems. In practice, F5+ terminates before reaching the degree cutoff, 
but it is not possible to test all systems, nor practical to determine a priori the precise degree of 
each Grobner basis. The question of whether F5, as presented in [16], terminates correctly on all 
systems, or even on all regular systems, remains an important open question. 

The following conjecture arises from an examination of Table [T] 

Conjecture 32. The F5 algorithm can terminate once all GB-critical pairs are rejected by the F5 
criteria. That is, it can terminate once d = (Ifr- 

Conjecture [35] is not a Corollary of Theorem [TJ] There, correctness follows only if all critical 
pairs are rejected by the algorithm: GB- and F5-critical pairs. Similarly, a proof of Conjecture |3"51 
would imply that we could drop altogether the check of Buchberger's criteria. 

If one could show that <i ma xGB < oIfr, Conjecture 1321 would follow immediately. However, such 
a proof is non-trivial, and lies beyond the scope of this paper. The conjecture may well be false 
even if we replace dp r by dp, although we have yet to encounter a counterexample. The difficulty 
lies in the possibility that Situation |2~T1 applies. 
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